Transport in the metallic regime of Mn doped III-V Semiconductors 
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The standard model of Mn doping in GaAs is subjected to a coherent potential approximation 
(CPA) treatment. Transport coefficients are evaluated within the linear response Kubo formalism. 
Both normal (NHE) and anomalous contributions (AHE) to the Hall effect are examined. We use 
a simple model density of states to describe the undoped valence band. The CPA bandstructure 
evolves into a spin split band caused by the p — d exchange scattering with Mn dopants. This gives 
rise to a strong magnetoresistance, which decreases sharply with temperature. The temperature (T) 
dependence of the resistance is due to spin disorder scattering (increasing with T), CPA bandstruc- 
ture renormalization and charged impurity scattering (decreasing with T). The calculated transport 
coefficients are discussed in relation to experiment, with a view of assessing the overall trends and 
deciding whether the model describes the right physics. This does indeed appear to be case, bearing 
in mind that the hopping limit needs to be treated separately, as it cannot be described within the 
band CPA. 

PACS numbers: 75.50.Pp, 72.10.-d, 72.25.Dc, 75.47.-m, 72.80.Ng, 71.70.Ej 



I. INTRODUCTION 



Magnets, which can be made using semiconductors by 
doping with magnetic impurities, are potentially very 
important new materials because they are expected to 
keep some of the useful properties of the host system. 
One important question is how much of the "doped 
semiconductor" bandstructure does the system still 
have? Is it a completely new alloy or does the system 
still behave like a doped semiconductor, which can, for 
example, be described using Kane theory, as taught in 
standard textbooks^. Many of these questions have been 
seriously investigated in extensive recent reviews2i2^£. 
The objective of this paper is to present a simple and 
mostly analytically tractable description of the mag- 
netism and transport properties of GaMnAs. We need 
a theory which captures the essential physics, can be 
parameterized and applied to magnetic nanostructures 
without having to carry out ab-initio calculations each 
time some parameter has changed. Previous theories 
of GaMnAs and related materials have been either 
computational or semi-analytical. But the treatment has 
often been split into three parts: magnetism, resistivity, 
and magnetotransport have been treated separately. 
Here we show that the magnetism and transport can be 
treated on the same footing using a one-band model. 
The one band model is of course not enough to explain 
the optical properties. But we show that magnetism, 
spin, and charged impurity scattering can be formulated 
in the same framework, using the coherent potential 
approximation (CPA) and a generalization thereof. In 
this paper we focus, however, only on the magnetism 
and spin scattering aspects and leave the computation of 
the charged impurity scattering and the full quantitative 
comparison to a later paper. One of the most difficult 
and controversial aspects for all conducting magnets is 
the origin of the Anomalous Hall effect (AHE). This is 



why we have devoted a section to reminding the reader 
of the basic definitions and results. The full history of 
the AHE is given in Ref. [H 

We focus on GaMnAs as a much-studied proto- 
type. It is now generally accepted that the observed 
ferromagnetic order of the localized spins in Mn-doped 
GaAs is due to Mn impurities acting as acceptor sites, 
which generate holes. These are antiferromagnetically 
coupled to the local 5/2 Mn spins, lowering their energy 
when they move in a sea of aligned moments. It is 
the configuration of lowest energy because spin-down 
band electrons move in attractive potentials of spin-up 
Mn and vice versa. Given that the spins are randomly 
distributed in space, the magnetic state is also the state 
of maximum "possible order". Thus, free and localized 
spins are intimately coupled. These materials exhibit 
strong negative magnetoresistance because aligning 
the spins reduces the disorder, and this lowers the 
resistance. The degree to which the magnetization of 
the spins affects the scattering process is dependent 
upon the degree to which the magnetic spin scattering is 
rate-determining for resistance. We shall, in this paper, 
only include the spin scattering process in order to first 
gain an intuitive understanding of the processes involved. 

The "hole mediated" magnetism point of view is 
not shared by everyone. Mahdavian and Zunger- argue 
that the mobile hole induced magnetism cannot explain 
the magnetism in high band gap, strongly-bound, hole 
materials such as Mn doped GaN. Their model, based 
on first principle supercell calculations, predicts that 
the Mn induced hole takes on a more d-like character as 
the host band gap increases, making the magnetism a 
" d — p" coupling property. 

The material of this paper is structured as follows: 
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we first present a short review of the phenomenology 
of the Hall effect in magnets. The Hamiltonian which 
describes the properties of Mn doped semiconductors is 
then introduced. We then discuss, in more detail, how 
to formulate transport in magnetically doped semicon- 
ductors. Following that, we recall how one can calculate 
the self-consistent CPA self-energy caused by the spin 
dependent term in the generally accepted Hamiltonian, 
for a one band system. Once the self-energy is known, we 
compute the longitudinal and transverse conductivities 
using the Kubo formulae. When the Fermi level is just 
above the mobility edge of the hole band, the Kubo 
formula is still valid and describes the strong-scattering 
random phase limit. The localized hopping regime is 
treated in another paper—. The CPA equations only need 
the density of states of the active free hole band as input. 
We will therefore consider the semicircular Hubbard 
band which reproduces the correct free hole band edges 
and is mathematically convenient for illustrating the 
effects of band induced magnetism. 



II. THE HALL EFFECT 

The experimentally measured Hall coefficient Rh is 
often written as 

Rh = Rn + Ra 

Ra = (a Pxx + bp 2 xx ) M/B z , (1 ' 

where a,b are constants and p xx is the resistivity. The 
first term Rn is the normal Hall coefficient and scales 
with resistivity in the usual way and the second Ra is the 
anomalous term, which in general can have two compo- 
nents, one linear and the other quadratic with resistivity 2 
and is proportional to the magnetization M. The general 
relation for the Hall coefficient is 

R H = Re W^)} 

B.^c^^,)}} 2 ' 

where a xy and a xx denote the transverse and normal 
conductivity respectively and B z is the magnetic field. 

Karplus and Luttinger— pointed out that the B- 
field involves the magnetic moment of the material via 
the internal magnetization M 

B z = Mo [Ht xt + (1 - N)M] = B e z xt + fi Q (l-N)M, (3) 

where N is the demagnetizing factor. Thus the mag- 
netization term is implicit in the normal contribution 
as a shift in the magnetic field. However this form 
is not normally sufficient to explain the much larger 
magnetization contribution observed in ferromagnets 8 . 
Throughout we shall, for simplicity, suppose a thin film 
with perpendicular-to-plane magnetic field and thus take 
N = 1, unless otherwise mentioned. 



III. MICROSCOPIC DESCRIPTION 
A. Hamiltonian 

In this section we formulate the basic model. The 
Hamiltonian for this Mn doped " alloy" , which most work- 
ers in the field accept as the correct description 2 -, is given 
by 

Htot = Hp + HLoc + Hd + Hpd + Hp-B+Hd-B + Hso- (4) 

In Eq. ((4]), the first term describes the free band holes. 
In the local Wannier or tight binding representation, it 
is 

Hp = ^ t tmn c ms c ns> (5) 
m,n,s 

with t mn denoting the overlap or jump terms from sites m 
to n, the c|j S) c ns are creation and annihilation operators 
for a carrier of spin s at site n, and where the indices 
include multiband transfers, if any. In the presence of a 
magnetic field, the overlap has a Peierls phase, so that 
we write it as (assuming thin film geometry, no internal 
field effect for a field perpendicular to plane) 

where ext stands for external, e is the electronic charge 
and R„ the position vector of site n. The second term 
in Eq. describes the diagonal (atomic orbital) energy 
that the valence band p-hole experiences when it is sitting 
on an impurity site which may or may not be magnetic 

Hl oc B ma c^ ra s c ms , (7) 

m,s 

where a stands for magnetic (M) or non-magnetic 
(NM). 

The third term is the direct exchange coupling be- 
tween the Mn d-localized spins, which we neglect here 
because when mobile carriers are present, they mediate 
the exchange coupling and in the case of a low concentra- 
tion, the Mn spins should be, on average, far away from 
each other. The fourth term is the antifcrromagnetic 
spin exchange coupling between the valence p-hole and 
the local d-Mn and ultimately the reason for magnetism 
in these materials ^ 10 ' 11 

Hpd = ~Y ^2 y^ C TOs°"^' " S TO<W) (8) 

m£{N s } s,s' 

where <r ss / is a vector containing the Pauli's matrices i.e. 
[a x ,a y ,c 2 ]. The indices, s and s' indicate which terms 
of the 2x2 matrix we are considering. Finally, S m is the 
Mn spin operator at site m. The fifth and sixth terms 
are the Zeeman energies of the holes and Mn spins in an 
external magnetic field along the z-axis, B z xt . 

H P -B = - yMB Yl B t% c lns<r Z s S Cms, (9) 
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Hd-B — -2/i B 2_j B t%Smi 
me{N B } 



(10) 



where g* is the effective g-factor. Finally we have the 
spin-orbit coupling 

H so = ^2^2(s,m\h so \n, s^c^Cns' , (11) 

m,n s,s' 



where 



4m 2 c 



(W(r) x p) 



(12) 



and where V(r) is the total potential acting at a point r 
and will include both normal crystal host sites and im- 
purity sites and p is the momentum operator. Finally, 
c is the speed of light. In the tight binding representa- 
tion, disorder is normally in the diagonal energies. For 
the spin-orbit term, disorder will enter the Hamiltonian 
through variations in the site potential. 



B. Spin-orbit coupling 

It is generally believed that the spin-orbit coupling is 
the cause of the anomalous Hall effect. This is one of 
the most fascinating and universal observations made in 
conducting magnets. In this section we briefly discuss 
the basic phenomenology of spin-orbit coupling. It is 
useful to recall the basic premises. 

When the Bloch wavefunction ansatz is inserted into the 
Hamiltonian, the spin-orbit coupling produces two new 
terms which enter the Hamiltonian for the periodic part 



of the wavefunction "it n k" 4 r nk(r) = Un^e 



Am 2 c 2 



(W(r)xp)-er + ft(W(r)xk)-er . (13) 



The first is the usual term and enters the band structure 
calculation. The second is called the Rashba terrn^ and 
in a crystal is, in general, less important than the first, 
because in the nucleus, momenta are larger than the 
lattice momenta k. 

The "Rashba term" is not normally treated in the 
computation of u n k. But in a lattice, it may be 
enhanced. In the Kane model, Chazalvieli 3 - and De 
Andrada e Silva et alJ^ have shown how to handle the 
spin-orbit terms in a lattice and the Rashba term in 
the presence of a triangular potential produced by a 
gate in an inversion layer. In particular, Ref. have 
shown how to renormalize the coefficient of the Rashba 
term and that this term acts when inversion symmetry 
is broken in an external field or triangular potential, 
for example. The origin of the enhancement of both 
these spin-orbit effects can perhaps be understood as 
similar to the origin of the change of the effective mass 
in the lattice. The spin-orbit energy can be enhanced 



by the presence of the lattice, but the enhancement can 
be very different from case to case, and needs to be 
examined in each case separately. For example, naively 
speaking, one may consider the spin-orbit force acting on 
a particle moving in a slowly varying potential of longer 
range than the lattice spacing as a spin current of an 
effective mass electron. One can think of the scattered 
particle as having an effective mass m* and generating 
a spin-orbit interaction which scales as instead 



of 



4 (rac)^ • When the particle is scattered from an atomic 



size impurity in a lattice, it is forced to come back many 
times and spends a longer time on the impurity than 
in free space. Leaving aside these intuitive pictures, 
one can, in any case, make the analysis quantitative 
using the Kane Hamiltonian, which takes the lattice 
modulation into account in a kind of renormalizcd 
perturbation theory and gives explicit results for the 
Unit part of the wavefunctionA. The Kane Hamiltonian 
can then be used to treat scattering from an impurity. 
Thus, for example, the matrix elements of the position 
operator (k|r|k') are very different for plane wave states 
and for the Kane solution s 8 ' 13 . 

In tight- binding (TB), the spin-orbit matrix ele- 
ments are calculated using atomic orbitals, with two 
terms: the intra- and the inter-atomic contributions. 
The magnitude of the intra-atomic term is just what it 
would be for the corresponding orbitals on the atoms 
in question. The inter atomic term is normally small 
and not included in tight-binding band structure cal- 
culations. When an electron jumps from site to site, it 
experiences a net magnetic field produced via the cross 
product of its velocity and the electric field gradient of 
the neighboring ions. This field then couples to its spin. 
Since it is related to the intersite momentum or velocity, 
it depends on the transfer rate t mn where the indices 
m,n include a band index m = m, 7. Remember that 
the velocity (x-direction) operator is given by 



. ^ (Rm,n) x t 



mn^ms^ns , 



(14) 



m,n,s 



where R m is the position vector of site m and 
Rm,n = R-m — R«- An enhancement of this spin-orbit 
energy, if any, has to be calculated by taking the expec- 
tation value of the spin-orbit term using the calculated 
Bloch states as one would for the dispersion relation Sk s 
and for the effective masses. We shall examine the spin 
orbit coupling in more detail below. 



IV. THE KUBO TRANSPORT EQUATIONS 

In order to compute the transport properties of the 
magnetically doped semiconductors in the linear response 
regime, we need to introduce the Kubo formulae. In this 
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section we show how to compute the transport coeffi- 
cients. 



A. General considerations 

The conductivity in linear response to an electric field 
is usually written a a 15 i 16 ' 17 



ihe 2 



lim 

<5->0 



E 

a,/3 



(a\v^\l3)((3\v^\a) f(e a )-f(ep) 
e a — ep + hu> + iS ep — e a 



(15) 

where the are the velocity operators in the respective 
direction ji. For a given Hamiltonian they are derived 
from the Heisenberg relation ihv^ — [x^, H] with x^ the 
position operator. The \a) and e a are the exact wave- 
functions and energy levels. The wavefunctions include 
any disorder and magnetic and spin-orbit couplings, /(e) 
is the Fermi function and u> is the frequency of the applied 
electric field; e is the electronic charge. This formula can 
now be written in the particular representation selected, 
i.g. tight-binding (TB) or Bloch states. Within the band- 
structure approach for semiconductors, one would sub- 
stitute the 8-band k • p wavefunctions and energies into 
Eq. (fl~5|) , include spin splitting through a Weiss field and 
spin-orbit coupling, and treat disorder as a lifetime con- 
tribution in the energy levels. In this picture the doped 
semiconductor retains the pure band structure features, 
apart from a complex lifetime shift. We shall use the one 
band CPA for the coupling of the holes to the magnetic 
impurities. For the calculation of the ordinary transport 
coefficients we may neglect the spin-orbit coupling. We 
obtain, for the longitudinal conductivity of the one band 
TB model for lu = and per spin s 



nil 



dE ) 



J deX(e) 



where 



(G s (E,s)) = 



lm{(G s (E,e))} 
(16) 



(17) 



and 



£ s — ^s[R] - i^s[I], 



(18) 



(19) 



Ti s (E) is the CPA self-energy for spin state s, to be calcu- 
lated self-consistently using the CPA condition and where 
the velocity and effective masses of the free electron band 
are given by 



1 de k 
h dk.. 



(20) 



1 d 2 e k 
ti 2 dkndk u ' 



(21) 



where M^ v is the effective mass tensor. The CPA 
self-energy contains the information on the net spin 
splitting produced by the magnetic impurities and the 
scattering time. It defines the effective medium which is 
to be discussed below. 

The Hall conductivity, given by the antisymmetric 
part of the transverse conductivity—^, is somewhat 
more complicated because the magnetic field has to be 
dealt with first. To first order in magnetic field, and 
with one TB band we find for the antisymmetric part 
(referred to by the index a) for oj = 15 i 17 ' 18 : 



2e 3 h 2 B 
3nfl 



where 



E 



dE 



dfjE) 
dE 



dsY{e) 



k L 



M. 



yy 



2VxVy 
M xy 



Im{(G s (E,e)) 
(22) 

5(e-e k ). (23) 



In Eq. ([22]) . the dispersion relation e k can be any one 
band structure (any choice of i^) where e k is defined 
through t® nn — jjj^k e lk Rm "£fc. As was shown by Mat- 
subara and Kaneyoshii^, for weak magnetic field, the 
Peierls phase can be factored out of the Green func- 
tion and this is why we only need the Green function 
(Eq. (fl"7|) ) that does not include the orbital effect of the 
field. The effect of the different phases (from the t mn and 
from the Green functions), after an expansion to first or- 
der in the field, are reflected in the topology dependent 
term Y(s). Furthermore, Eq. (|22|) does not contain the 
spin-orbit scattering contribution. It can be included as 
an extrinsic effect. In the spirit of the skew scattering 
Born approximatio n 19 i 20 i 21 and in the weak scattering 
limit, we can, at the end of the calculation, introduce 
an extra internal magnetic field so that the total field 
entering Eq. (f2"2"| is 



B = B Z 



m 
e 



(24) 



where B z includes the internal magnetization. However, 
from Eq. ((3]), for a field perpendicular to plane, N = 1. 
This result, Eq. (|24| . is written in the notation of 
Ballentinei 9 -. The effective spin-orbit field depends on 
the thermally averaged spin polarization (cr z ), with 

its sign to be discussed in Section IIVB1 and {^7-^) is 

the so-called skew scattering rate. The second term of 
Eq. (|24[) involves, in addition to the sign of the carrier, 
the sign of the spin orientation. 

Engel et al. 21 claim that the impurity spin-orbit 
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coupling in the lattice environment can be "6 orders or 
magnitude" larger than in vacuum. The large enhance- 
ment can mean that, in some cases, the skew scattering 
dominates where there should exist a substantial in- 
trinsic contribution as wel l 21 i 23 . This is why developing 
a way to calculate the extrinsic contributions to the 
anomalous Hall effect, as we do here, is important and 
useful. 



B. The Anomalous Hall Effect (AHE) 



A„(r) 



4m 2 c 2 



W„(r-R„) 



r R n 



(26) 



where A n (r) is related to the spin-orbit coupling 
strength. Here Pi nter i s the usual zero-field inter-atomic 
momentum operator, lj is the atomic orbital operator 
and is often quenched so that its expectation value is 
zero, but nondiagonal same site inter-orbital matrix 
elements can be very important. 



Let us now show how one can derive such an anoma- 
lous Hall effect, in our formalism (TB approximation), 
from the existence of a spin-orbit coupling and how one 
can arrive at the concept of an effective magnetic field, 
as given in Eq. ([24)1 . 

In tight-binding, the spin-orbit coupling is given 
by Eq. and V(r) is the total potential experienced 
by the charge at point r, with p the momentum operator. 
The spin-orbit coupling can be separated into an " intra" 
and an "inter" atomic contribution by splitting the 
momentum operator as p = p a + Pi„ ter ■ For spherically 
symmetric potentials the two terms are given by 

Vso =^Ai(r)lj • er + ^ j(r - R„)A„(r) x p inter j ■ a, 

i n 

(25) 



The second contribution in Eq. ([25]) is due to the 
inter-atomic motion, and must be evaluated using 
Pmter = mv where v is given by Eq. (fT4|) and involves 
the intersite transfer energy t and position operator 



I" — ^ ] 'R-mci ns c ms- (27) 



Here we have a spin-orbit coupling only because the par- 
ticle can jump to another orbital. When substituting the 
site representation in the second term of Eq. (|25p . and 
assuming spherical symmetric potentials, either impurity 
with effective local charge eZ n = eZi mp or host eZh we 
have 



(i\V s ^\j) = 



Amc 2 



E 

n,l 



(4tt£ ) |r - R„ 



r(r-Rn) 



(28) 



r 



where i is the imaginary number while (i\ is the orbital 
at site i. The analysis of this term is not trivial, and 
depends on the lattice and potential distribution in ques- 
tion. If one has decided that the extrinsic contributions 



are dominant, then one only needs to sum around the 
impurities. In general, the simplest approach is to take 
only the largest terms in the sum, i.e. the diagonal terms 
I = i. In this case we get 



J 



{i\V s i ?\j)=i(i\^\j)t ij , 
h 



i(i\®<r\j) = 



Amc 2 




(r - R„) 



1 H XR - 



cr, 



(29) 



It is convenient to treat this effect as a spin dependent can write i&ij = e l ® ij — 1 and H, p + H so as 
phase in the transfer, in analogy to the Peierls phase, to 

first order. Indeed, to first order e x rs 1 + x so that we „_ _ , 

i,j cr.cr' 



(30) 
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Then we also note that the R„ 



R, and R M 



R, 



terms in Eq. (j2T)|) are zero after the vector product, 
leaving the gradient field contributions due to the 
nearest neighbors as the largest of the remaining 
terms in the n-sum. A similar result has been derived 
in the context of the Rashba coupling by Damker et ali^I. 



With Eq. 
element as 



we can rewrite the spin-orbit matrix 



h s ° = i—t l3 ^ (Hj x R m ) • B n , s0 , (31) 

where Rj„ = R; — R„ and where we have introduced a 
spin-orbit magnetic field defined by 



B„, so (i -» j) = 



2mc 2 



Z,, 



47T£ |R 4 - RJ 



(32) 



In summary, let us write in this approximation, a 
new total phase for Eq. ^ 



<Pij 



2^ I — B z • (Rj x Rj) + (R-j x R-in) • B„ jSO 



We thus see what the spin-orbit coupling does to the 
TB Hamiltonian. It is, in terms of energy, a small ef- 
fect. Its effect on the band structure and magnetism will 
be neglected. In order to obtain Eq. (|22[) with B de- 
fined as Eq. ([24]) one as to decouple the trace over the 
spin polarization from the remaining expression and this 
automatically makes (a z ) the overall mobile spin polar- 
ization. In reality, the skew-scatterring is due to carriers 
near the Fermi level and therefore (<r z ) should be the av- 
erage weighted at the Fermi level. We will return to this 
point when we attempt to fit Eq. (j2"2"|) and Eq. (f2~4"|) to 
experiment in Section [VI Bl 



C. Transport coefficients for a simple cubic band 

To perform a calculation, one needs to specify a lat- 
tice topology (i^n)- The simplest one is a simple cu- 
bic lattice with nearest neighbor hopping. The disper- 
sion relation is given, for a d-dimensional system, by 
£k = — 2t 53 a =i cos(fc Q a) where a is the lattice constant 
of the simple cubic. With this particular form, we can 
show-2£ that the longitudinal and Hall conductivities at 

and 



zero frequency can be written, using Eqs. H 
with no approximations, as 



(cTxx) 



d7rfot d - 2 



E 



df(E) \ 
dE J 



111! 



{(Gs(E,e))} 



-zDo(z)dzdedE 



(34) 



4e 3 B 



37rd(d - l)a d ~ 4 h 2 



E 



df(E) \ 
dE J 



Im 



{(Gs(E,e))} 



-zD()(z)dzdedE 



(35) 



where (G S (E, e)) is given by Eq. (fTTj) . d is the number of 
dimensions, a stands once again for antisymmetric and 
D (z) is the density of states of the pure crystal. In 3D, 
Dq(z) can be found knowing that, for the above for 
a d-dimcnsional cubic lattice, the Green function is^ 



°o( E ) = Al / #*K(aO, (36) 



where K(x) is the complete elliptic integral of the first 
kind, with x = E+lST 2 tC osW 

Now we can return to the calculation of the new 
bands generated by the hole-spin coupling and for this 
we use the CPA. Later, we use Eq. (|24l) to estimate the 
magnitude of the skew scattering Hall effect. 



V. THE CPA EQUATIONS FOR A ONE-BAND 
MODEL WITH COUPLING TO LOCAL SPINS 

In this section we show how to compute the new energy 
bands in the presence of a high concentration of dopants, 
magnetic and nonmagnetic. The magnetic coupling 
between the Mn spin and charge and the holes in the 
valence band creates a new valence band-structure. This 
new Mn-induced band will be magnetic at low enough 
temperatures. Since the Mn are randomly distributed, 
we cannot use Bloch's theorem. One way forward is to 
use the powerful self-consistent single site approximation 
known as the CPA. The CPA self-energy 30 T,(E) is 
determined by the condition that if at a single site, the 
effective medium is replaced by the true medium, then 
the configurationally averaged t-matrix produced by 
scattering from the difference between the true medium 
and effective medium potentials must vanish 30 . 
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In the present case we have localized Mn spins of 
magnitude S — 5/2 so that there are 6 possible states 
S z = —5/2,..., 5/2 and thus 7 possible sites in the 
system including the non-magnetic ones. We have for 
a "spin up" carrier, the normal non magnetic sites 
with concentration (1 — x) with local site energy Enm 
and magnetic field interaction — ^-/Ub-BJ 3 ^) and the 
magnetic sites with concentration x with local site 
energy Em and magnetic coupling ^rS^ — ^-(IbB^^. 
But there is also the possibility of a spin flip of both the 
carrier and the impurity with interaction -§^5+. For 
a "spin down" carrier we have Enm at a nonmagnetic 
site and magnetic field interaction ^[IbB!^. For a 
magnetic site we still have Em as the local site energy 
but the magnetic coupling is now — ^pS^, + ^tUbB^^ 

and the spin-flip interaction is given by ^-S^. The 
concentration of spin carrying impurities is defined by 
x and is where Nq is the number of sites with an 
impurity while Nl is the total number of sites in the 
system. The CPA conditions are given b y 28 ' 31 



(l-x)(C%) Th + x(t% n ) Th = 
(l-*)<0™ + *<<U>™=°> 
and the thermal average of an operator is 



(37) 



{O m {S z )) = 



Tr (O m {S z )e 



Tr(e 



-0H S 



t3h m (S*)S* 



(38) 



In order to find the spin dependent t matrix operator 
for the CPA, it is useful to rewrite the Hamiltonian as a 
2x2 matrix in spin space. The t-\-\ and t±± operators are 
the diagonals element of the i-matrix. The expressions 
for the t-matrix are essentially the same as the ones ob- 
tained previously by Takahashi and Mitsui 3 ™. The minor 
difference is that in the present case, we include the mag- 
netic field, so that the Zee-man term enters the t operator. 
The inclusion of the Zeeman term is straightforward: one 
has to add, when N = 1, ^-urB^* with the appropri- 
ate sign in the expressions of Ref. [U To calculate the 
CPA i-matrices, we need the effective medium Green's 
function which is no longer dependent on the site m. 



G S (E) = (G mms (E)) = 



1 



N ^ E- e ks 

k 



E.(.E) 



(39) 



In the previous equation, Gij S (E) refers to the Green's 
function of a carrier of spin s between sites i and j for 
one particular configuration of disorder. The average is 
over all possible realization of disorder. 



The local density of states is 

D S (E) = -hhn.[{G mms (£?)>}, 



(40) 



and the hole spin concentration with spin s at site m can 
be written 



(p ms ) = / dEf(E)D s (E). 



(41) 



Assuming mean field theory for the energy entering the 
Boltzmann factor, the probability that the local spin has 
a value S z is 



P{S Z 



-l3h m S* 



-t3h m S* 



where 



ext 
z.m ■ 



(42) 



(43) 



To compute the CPA self-energy T, S (E), we need to spec- 
ify the lattice dispersion of the hole band. In effect, it 
suffices to specify the bare pure lattice density of states 
(Dq(E)) since the real and imaginary parts are simply 
related. The Fermi level is determined by the condition 
that the total hole concentration p is known and fixed 
relative to the total impurity concentration with x as 
maximum value. In general, there will be fewer holes 
than dopants, but the number is not known and remains 
a fit parameter, so we have 



P = y dEf(E)[D 1 (E)+D l (E)}. 



(44) 



This now allows us to compute the CPA self-energy self- 
consistently and then to determine the new density of 
states via Eq. (|4"0")) and the local spin polarization and mo- 
bile spin concentration via Eq. (|38|) and Eq. (|4T|) . Thus, 
we can determine the magnetization and the transport 
coefficients, the resistivity and the conductivity as a func- 
tion of B via Eq. (f3"4"|) , and the magnetoresistance as the 
relative change with B. Finally, the normal and anoma- 
lous Hall Rh can be obtained via Eq. 



VI. APPLICATIONS OF THE CPA 

For numerical calculations, to simplify the analysis and 
for proof of principle, instead of the Dq{E) found with 
Eq. ((551) j we use t ne Hubbard function, given by 

nr (E) = sta^, (45) 

irW z 

where W is half the bandwidth and Q is the Heaviside 
step function. Eq. ([45)) has the same band edge behavior 
as the D (E) calculated with Eq. (|3"rJ)) . With the simple 
form of Eq. (|45p. the first integration in Eq. (f3"4")) and 
Eq. (|3"5)) is analytically tractable. 
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A. dc conduction and magnetoresistivity 

Here we calculate the dc conductivity {(J X x) from 
Eq. (j3~4"|) and the magnetoresistance. The calculations 
have been done assuming N = 1 (thin him). The calcu- 
lations for the magnetoresistance have been performed 
without taking into account the spin-orbit interactions. 

Figure [1] and Fig. [2] are plots of the CPA density 
of states for two values of p and J p d assuming that 
the scattering is solely due to the spin potential, that 
is when Em = -Ejvm = 0. The diagrams shows the 
evolution of the spin dependent density of states as a 
function of temperature for four different temperatures 
in each figure. The temperature is measured in units 
of the bandwidth W. In Fig. |l(d)| and |2(d)| one can no 
longer see the spin-splitting. The position of the Fermi 
level is also shown. Note that all energy parameters 
are normalized by W (E' = E/W). From Fig. Q] and 
Fig. [U one can see that when the magnetic coupling J v d 
is very small, the split band disappears (see Ref. [28| for 
details at low J p d), whereas in the opposite limit, the 
spin bands split off and a pseudo gap-appears. 

Figure. [3] shows the resistivity in zero B field as a 
function of temperature with p as a parameter. Results 
for two values of J p d are presented. The re sults for 
lower J p d than 0.35VF* are similar to Fig. 3(a)[ !g . When 
the impurity bands completely splits from the valence 
band ( for example for J p d = 0.5W ), the qualitative 
behavior changes, as there are now two well defined 
bands, above and below the Fermi level (impurity and 
valence) separated by a gap (see Fig. [3] above and 
Ref. [28| ). At high hole concentration (p > 0.8x), the 
system is metallic and resistance increases with spin 
disorder, at first rapidly, and then decreases again above 
T c . This is intuitively to be expected because at first, 
as temperature increases, the disorder increases, and 
then in the paramagnetic phase, thermal broadening 
overcomes the potential scattering disorder and the 
lifetime averages out. We should remember that in the 
simplest Boltzmann approach the conductivity is given 
by 



df(E)\ . 

-be- r* 



(E) Ts (E)dE (46) 



and at low temperatures, depending on the product 
of the scattering time and the density of states at the 
Fermi level. In an alloy, either quantity can change with 
B and T and determine the conductivity. Above T c the 
spin splitting disappears and the quantities involved 
are, in the absence of charged impurity scattering, 
only weak functions of temperature. At very low hole 
concentration, the Fermi level is in a region of small 
density of states where we expect localization. However, 
CPA is a mean field method and does not produce 
localization. Even though the resistance increases with 
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FIG. 1: Density of states (DOS) for different temperatures, 
(a) k B T = 0, (b) k B T = 5xlO" 3 W, (c) k B T = 7.99xlO" 3 W 
and (d) k B T = 8.1xlO~ 3 VF, when x = 0.053, E NM = E M = 
0, fi B B = 0, Jpd = 0AW and p = 0.8a;. The full line curve 
is the DOS for the spin down carrier, the dashed one is the 
DOS for the spin up and the vertical line shows the Fermi 
level. One can see that the density of states is, relative to a 
nonmagnetic system, a strong function of temperature. 
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FIG. 2: Density of states (DOS) for different temperatures, 
(a) k B T = 0, (b) k B T = 2.5xl(T 3 W, (c) k B T = 4.5xlO" 3 W 
and (d) k B T = 5.4xlO _3 W, when x = 0.053, E NM = E M = 
0, /j, b B = 0, J p d = 0.5W and p = 0.3a;. The solid line curve 
is the DOS for the spin down carrier, the dashed one is the 
DOS for the spin up and the vertical line shows the Fermi 
level. One can see that the density of states is, relative to a 
nonmagnetic system, a strong function of temperature. 
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FIG. 3: Resistivity as a function of temperature for two J p d, 
(a) Jpd = 0.35W and (b) J pd = 0AW for an impurity con- 
centration x = 0.053, Emm ~ Em = and \i B B = 0, for two 
values of J p d , with p as a parameter. 



decreasing density of states at the Fermi level, the 
conduction process continues to be band conduction 
albeit with short relaxation time. Even if the mean 
free path reaches the lattice spacing (random phase 
limit), the conductivity is still far greater than hopping 
conductivity between localized levels^. The effect 
of temperature on resistivity above T c is not too sig- 
nificant in the non-localized intermediate cases. The 
temperature dependence of the resistance is, in general, 
a complex interplay of density of states, velocity and 
relaxation time (imaginary part of the self-energy). The 
temperature induced lowering of the charged impurity 
screening length (see Eq. (|49[) to Eq. (|5Tj) ) is not included 
here. In this range of exchange coupling, J p d does not 
seem to strongly influence the structure of the resistivity 
versus T curves but does change their magnitude. In 
Fig. [31 for p = O.lx, a low carrier concentration, one 
may see that the resistance increases with J p d- In 
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FIG. 4: Resistivity as a function of impurity concentration 
x, for varying degrees of hole doping. The others parameters 
are J pd = 0.35W, E NM = E M = 0, k B T = 3.5xl0 _3 W and 
HbB = 0. 



FIG. 5: Resistivity as a function of temperature for various 
magnetic fields. Others parameters are x = 0.053, Enm = 
Em = 0, J p d = 0-4W and p = 0.8a;. The solid line is the 
zero field result, the dashed line is for hbB = lxlO -5 !^, the 
dotted line is for \ibB = lxl0 _4 VF and the dashed-dotted 
line is for fj, B B = 3xl0~ 4 W/. 



contrast, in the high hole concentration limit p = x, 
there is only a relatively weak variation of resistance 
behavior with J p d- The complex but regular behavior 
of resistance with temperature is a manifestation of 
Eq. (|46|) . reflecting the different elements which deter- 
mine the value of resistance for a given set of parameters. 

Figure [4] shows the resistivity at one temperature, 
as a function of impurity concentration for varying de- 
grees of hole doping. The results of Fig. [3] and Fig. 2] are 
based on the spin scattering model with no additional 
potential scattering terms, and no elect ron-phonon 
interaction. A fit to experimental data must take these 
other mechanisms into account as well. Thus, the 
complete CPA self-energy should include other sources 
of random potential scattering (nonzero values of Enm, 
Em), in particular, the charged impurity scattering 
terms discussed below. We should also add, when 
necessary, the electron- phonon self-energy T, ep (E). The 
imaginary part of the charged impurity self-energy will 
contribute another source of temperature dependent 
lifetime broadening, and the real part will enter the 
density of states. 

Figure O shows the resistivity as a function of tem- 
perature for various magnetic fields, for one value of x, 
p, and Jpd, and again with Em — Enm — 0. Apart 
from very low and very high temperature, increasing 
magnetic field decreases the resistivity, by favoring the 
alignment of the impurity spins. The magnetic field also 
eliminates the sharp metal-insulator transition around 
T c , replacing it with a smooth transition. 

Figure O shows the relative magnetoresistivity de- 



fined by 



MR ^P*AB ext )-p x AB ext = 0) 



pxx{B e 



0) 



(47) 



as a function of B, for various value of T. The overall 
trend in this exclusively spin scattering model is strong 
negative magnetoresistivity, as one would expect, since 
increasing Mn spin alignment reduces disorder, and thus 
reduces the scattering lifetime, and no other source of 
scattering are considered. However, lifetime is not the 
only quantity entering the conductivity. As shown in 
Eq. (j46|) . the density of states D(Ep) at Ef also plays 
an important role. There are also regimes of positive 
magnetoresistivity. From Fig. 6(a)| we see that it is also 
possible for the resistance to increase with B at low T 
when the hole density is very low. This is probably be- 
cause in this limit, the density of states at the Fermi level 
decreases with magnetic field and this effect is stronger 
than the concomitant increase of the carrier lifetime due 
to the suppression of spin disorder with B. But we also 
know that for low concentrations, when the Fermi level is 
at the band edge or in the region of localized states, other 
changes arise which are not due to spin disorder scatter- 
ing and which require another approach which is based 
on localization. In the hopping regime, not describable 
by CPA, the magnetic field can, for example, squeeze 
the localized wavefunctions and increases resistance by 
reducing overlap. But it also shifts the energy and the 
mobility edge such as to reduce resistance^. When the 
calculated resistance is high and when the Fermi level is 
in a region of small density of states < 10 19 /cm 3 eV, we 
should therefore not trust the CPA, which gives a good 
description only of the metallic regime. 



11 





2 f 




n 




-2 




-4 




-6 






tr 


-8 






-10 - 




-12 




-14 




-16 



kT= 3.64x10 4 W 

. B . . 




k T=0.015W 
B 




k B T= 1.31x10 ~ 3 W 




k B T= 1.80x10 " 3 W . 







0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 
H B B/W (10" 4 ) 



k„T = 3.64x10 W 




|i B B/W (10" 

(b) 



FIG. 6: Relative magnetoresistivity for two carrier concentra- 
tion, (a) p = 0.1a; and (b) p = 0.8a:, for impurity concentra- 
tion x = 0.053, Enm = E M = and J pd = 0.35W, with T as 
a parameter. 



B. Hall effect 

By definition, the Hall resistance is given by 
-Rejo-yz j 



Pyx 



Rcjcr^j - Rejcr^ jRcjcr^ j 



(48) 



Equation ()48|) being a general definition, the normal and 
the anomalous spin-orbit generated Hall resistance are 
included. 

In the CPA, the sign of the normal Hall effect is 
not determined by a simple relationship, but depends 
on the dispersion of the lattice and the behavior of the 
real part of the Green's function at the Fermi leve l 15 ' 28 . 
A full analytical analysis of the CPA Hall sign is beyond 
the scope of this article but we can state that there is no 



simple rule. The closest to one is that Rh is electron-like 
when the density of states increases with energy and 
hole-like when it decreases. In the approximation of 
skew scattering, the sign of the anomalous effect follows 
the sign of the normal effect as long as the carriers at 
the Fermi level are polarized in the same direction as 
the total magnetization. This is true here, as can be 
seen from Fig. [T] The majority spin band at the Fermi 
level has the same sign as the overall magnetization, 
which is dominated by the localized spins. Thus, in the 
phase addition approximation of Eq. (|33[) . in order to 
get the right sign, it is essential to keep the correlation 
between spin and energy. The decoupling of the spin 
magnetization out of the Kubo formula can give rise to 
the wrong sign. Here, the overall polarization, found 
using Eq. (|4"Tj) . is opposite to the direction of the field as 
it should be. We correct for this effect by assuming that 
(a z ) is in the same direction as the overall magnetization. 

Figure [7] shows, for a selected class of parameters, 
the overall behavior of the Hall resistivity in the skew 
spin scattering model. In order to explain the exper- 
imental data we need — ttt = 0.2, which means that 

T s At ' 

the skew scattering energy has to be 0.4 times the 
tight binding overlap. This is a very (too) strongly 
enhanced skew scattering rate. It suggests that the 
correct interpretation of the AHE in GaMnAs is most 
likely the intrinsic mechanism proposed by Sinova et 
alA In this work, the AHE is due to the intrinsic 
spin-orbit field, and relies on the multiple band nature 
of this class of semiconductors. The usual simple 
nearest neighbor tight-binding model only gives a skew 
scattering contribution. The universality and order of 
magnitude of the AHE in ferromagnets suggests that 
the intrinsic process dominates in most cases. 



At high temperature, the 
The normal Hall conduction, 

covered (see the ksT — 0.015VF results Fig. 7(a) and 



magnetism disappears, 
linear in B field, is re- 



Fig. |7(b)[ ). Note that for simplicity of notation, we refer 
to the applied field in the figure as B, even if it was 
defined otherwise previously. 



VII. EXPERIMENTAL RELEVANCE 

Figure [8] shows results of the CPA calculations for the 
three transport parameters: resistance, magnetoresis- 
tance and Hall effect in a range of parameters for which 
a behavior close to the ones observed experimentally by 
Ruzmetov^ and Ohno et al. 35 is observed. We have 
added a constant contribution to the resistivity so that 
the relative change in resistivity, such as discussed in 
connection with Fig. between T c and T — is of the 
same order of magnitude as observed in experiment. 
The constant term is chosen to be of the same order of 
magnitude as the spin scattering rate. The qualitative 
temperature structure is satisfactory in the 'metallic 
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FIG. 7: Hall resistivity as a function of the applied magnetic 
field for two carrier concentration, (a) p = O.lx and (b) p = 
O.81, for an impurity fraction x — 0.053, Enm = EM — 
and J p d = 0.35W, with T as a parameter 



regime' of the material, but because we have neglected 
the charged impurity scattering, we are underestimating 
the temperature drop of the resistivity with tempera- 
ture, at higher temperatures. The latter is due to a 
reduction in screening length when the density of states 
at the Fermi level increases, (see Eq. (|5Tj) ). Figure [5] 
also illustrates the strong temperature dependence of 
the magnetoresistance, which reaches 40 to 50% at 
k B T ~ 0.0057W (this is - 60X for W = leV), and 
how it is intimately connected to the magnetic order as 
observed experimentally. For an even clearer picture of 
the connection see Fig. [5] The reader should note that 
Van Esch et al^ have observed magnetoresistances of 
~ 500% at T = AK and 50% at T = 20K in their more 
resistive samples with p ~ lf2/cm. In this sample the 
magnetization dropped by ~ 100% between AK and 20JT. 

Hwang and Das Sarma^i and Lopez-Sancho and 
Brey2^ quite rightly pointed out that when charged 
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FIG. 8: Transport coefficients which can be directly compared 
to experiment, (a) p = po + p(T), (b) MR(B) and (c) ph(B), 
for E NM = E M = 0, J pd = 0.35IT and p = 0.8a;. Two 
possible constant contributions are shown in (a), the solid 
line curve has po = lOmQcm and the dashed line curve has 
po = 3.6mf2cm 



13 



impurities are present, they will normally dominate 
the scattering rate and the temperature dependence 
of resistance in a wide range of parameters. They 
then decided to use two different conductivity models 
for the transport in Mn doped GaAs, one to explain 
the temperature dependence of the resistance which 
emphasizes charged impurity scattering, and one for the 
magnetoresistance which emphasizes spin scattering. 
We believe this is not necessary. It is enough to add the 
configurationally averaged charged impurity self-energy 
to the basis band dispersion £ ks ~~ * £ks + Sk s , and then 
to neglect, in the simplest limit, the real part of the 
impurity scattering self-energy £ ks • Thus we replace the 
CPA self-energy, S s , which appears in Eq. (fTI)]) . by 

Im{E s (£)} +Im{E ka (£)} - Im{s s (£)} + 

r, J dk' |V k _ k f (1 - cos0 k _ k ')<5(£ - e k / 

(49) 

where rj = jf-Ni, with N t the impurity concentration, 
and 



Vx 



k-k' 



AttZ- e 2 



1 



(50) 



££ |k-k'| l + (g sc /|k-k'|) 
In Eq. (|50|) . Zi mp is the impurity charge number, and 

ql = 4rr e 2 £ J dED s (E) (-^) (51) 

s 

is the inverse screening length. Since the results of Sec- 
tion EH show that D S (E) at the Fermi level change with 
temperature and field, the screening length will change 
as well. Note that usually, in Boltzmann transport, one 
does not take into account the effect of the real part of 
the self-energy (the imaginary part gives the scattering) 
on the actual band structure because the concentration 
of dopants is low. Here the concentration is not low. The 
real part should in principle, modify the band gap and 
effective masses and can be included in our formalism. 

We agree with Hwang and Das Sarma^ 7 - that the 
complete problem in parameter space is enormously 
complex. In effect, we should want the impurity scat- 
tering to only represent an additional self-consistent, 
density of states dependent, lifetime process. A fully 
self-consistent Coulomb potential/spin multiband CPA 
is far too complex, and of little value. In addition, 
the CPA band theory cannot account for the hopping 
regime, as pointed out in the previous section. 



VIII. THE PROBLEM OF THE INTRINSIC 
AND EXTRINSIC HALL EFFECT 

In the Kane-Luttinger theory of Jungwirth et al£, the 
Mn dopants only cause a lifetime broadening and no 



change to the Kane-Bloch bandstructure. In this formal- 
ism, the intrinsic spin-orbit interaction is Bloch invariant 
and only causes band mixing. All disorder is treated only 
as a lifetime effect. This is not so in the CPA where the 
dopant scattering causes a substantial renormalization of 
the bandstructure via the real part of the self-energy and 
indeed gives an explanation of the magnetism. The spin- 
orbit coupling produces new terms in the Hamiltonian 
which allow the electrons jumping or tunneling from one 
site to the next site to experience the electric field of the 
neighboring atoms. The 3-site processes we invoke are 
not normally included in the tight-binding modeling of 
semiconductors. They represent only small modifications 
of the bandstructure, much smaller than the on-site spin- 
orbit admixture. The anomalous Hall mechanism in one 
band tight-binding is basically one of " skew scattering" , 
except that the jump is always from orbit to orbit, and 
the effective magnetic field can only come from another 
3 rd site. We have seen that the one band TB description 
is inadequate to describe the AHE in DMS. On the other 
hand, the multiple band Kane and/or Kohn-Lutting 
approach seems to work well. It gives the right order of 
magnitude with an intrinsic AHE. This suggests that a 
TB modeling of the AHE must include the multiple band 
aspect and at least a second nearest neighbor overlap. 
Since, for example, a five percent GaMnAs alloy cannot 
have a true Bloch band structure, in the orbital descrip- 
tion, it is the short range many orbital aspect which must 
give rise to the intrinsic AHE. 



IX. CONCLUSION 

We have presented a CPA theory which could explain 
the magnetism in Mn doped semiconductors in the 
metallic regime and allows one to calculate the transport 
parameters. The numerical evaluation of resistance, 
magnetoresistance and Hall coefficient, normal and 
anomalous, show that the overall trends observed 
experimentally are reproduced by the CPA. The power 
of the method lies in its simplicity. All the information 
is contained in a spin and energy dependent self-energy 
Ti s (E). In good approximation, we may add the effect 
of the charged impurity and electron-phonon lifetime 
corrections to the calculated CPA self-energy. Of the 
two, the first is the most important addition because it 
explains why the experimental resistance decreases at 
high temperatures (in general more strongly than shown 
in Fig. 8(a)), when the ferromagnetism has reduced 
and spin-disorder scattering should 
The "small" resistance 



to paramagnetism, 
have reached its highest value, 
drop at high temperature seen in Fig. 8(a) is due to CPA 
bandstructure renormalization. The CPA and i-matrix 
methods can be extended to treat magnetic clusters and 
evaluate effective localized spin-spin coupling mediated 
by the band. 

We have also demonstrated how to include charged 
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impurity scattering within the same formalism. The 
central aim of this paper was to achieve an understand- 
ing of the important mechanism which determines the 
conductivity behavior. For this purpose, we have used 
a one-band approach which is adequate to understand 
the conductivity. For optical properties and the intrinsic 
AHE, the many band aspects are essential. 

The AHE has been modeled as an extrinsic effect 
in the framework of skew scattering. We did this using 
the tight-binding language. The intrinsic AHE, as 
invoked by Jungwirth et al£, cannot be derived using a 
nearest neighbor TB formalism, and it is not clear how 
to recover it in the TB formalism. The problematic of 
the intrinsic AHE and the TB model is an interesting 



one. It should be re-examined in detail because so far, 
the nearest neighbor tight-binding methods have proved 
useful as bandstructure descriptions of semiconductors. 
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